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A b s tra ct 

To assess the state of the art in ship 
viscous flow computation a Workshop was 
organized in 1990 by three organizations: 
SSPA Maritime Consulting AB, Chalmers Uni- 
versity of Technology and Iowa Institute of 
Hydraulic Research. Two test cases were 
specified by the organizers and sent out to 
all interested research groups, which were 
asked to submit results in a prescribed 
format. In Septemer 1990 a meeting was held 
at Chalmers University of Technology. All 
results had then been collected and presen- 
ted in a common format, and the theories 
behind the methods compiled in a table 
based on responses to a questionnaire sent 
our earlier. During the meeting, each 
research group was first given the oppor- 
tunity to briefly introduce their method 
and results. Thereafter, a considerable 
time was spent on general discussions on 
the performance of the different methods 
considering the differences in the under- 
lying theories. Specific items that were 
adressed were grid generation, governing 
equations, boundary conditions, turbulence 
modelling and numerical method. Practical 
aspects on the results, for instance from 
the point of view of propeller design, were 
also discussed. The Workshop Proceedings 
contain a description of the participating 
methods and the results of both test cases. 
In the present paper a summary of the Work- 
shop and its results is presented. 

Introduction 

Although viscosity is present in the 
entire flow around a ship it has a signi- 
ficant effect only in the boundary layer 
around the ship and the wake behind it. The 
present paper deals with the flow in these 
two regions, which are limited in size, but 
very important from at least two points of 
view. Frictional forces within the boundary 
layer give rise to a viscous resistance of 
the ship, in most cases the dominant resis- 
tance component, and the velocity distribu- 
tion in the near wake determines the design 
and performance of the propeller. 

Despite its obvious importance the 
first serious attempts to compute the 
viscous flow were made relatively recently, 
about twenty years ago. This is in contrast 
to the long term research in the inviscid 
flow area, where wave resistance research 
has been under way during the entire twen- 
tieth century. The reason for the diffe- 


rence is that the complicated viscous flow 
equations, i.e. the Navier-Stokes equations 
and approximations thereof, are less ame- 
nable to analytical treatment than the in- 
viscid equations, normally based on poten- 
tial flow theory. Therefore it was not un- 
til computers had become powerful enough to 
handle three-dimensional boundary layer 
theory numerically that research on viscous 
flow computations was started. 

During the 1970 f s a number of methods 
for predicting ship boundary layers were 
developed, and in 1980 it was considered 
appropriate to assess the state of the art 
in this area. To accomplish this, an inter- 
national Workshop on ship boundary layers 
was organized by SSPA, the Swedish State 
Shipbuilding Experimental Tank, in coopera- 
tion with the International Towing Tank 
Conference, ITTC 1 . The purpose was to bring 
together specialists on ship boundary layer 
calculations from all over the world, and 
to let them apply their techniques and 
methods to two test cases, specified in 
detail by the organizers. In June of 1980 a 
meeting was held in Gothenburg. The results 
of 17 methods had then been collected and 
presented in a uniform format. During the 
meeting the various components of the 
methods were discussed in the light of the 
results produced. The general finding was 
that most methods were able to predict the 
thin boundary layer over the major part of 
the hull with an accuracy sufficient for 
engineering purposes, while all of the 
methods failed completely in predicting the 
flow near the stern and in the wake. 

In the 1980 r s development accelerated, 
and the focus was changed from the thin 
boundary layer to the stern/wake flow. As 
evidenced by the 1980 Workshop a new class 
of methods with less restrictive approxi- 
mations was required. The interest of 
researchers was soon focused on the Rey- 
nolds Averaged Navier-Stokes (RANS) equa- 
tions, and a number of such methods was 
proposed during the 1980* s. Towards the end 
of the decade it was considered timely to 
organize a second Workshop to investigate 
the progress made. This task was undertaken 
by three organizations: SSPA Maritime Con- 
sulting and Chalmers University of Techno- 
logy (CTH) in Sweden and Iowa Institute of 
Hydraulic Research (IIHR) in the USA. 

As in the first Workshop, the purpose 
of the new event was twofold: 


o to assess the state of the art in ship 

viscous flow calculations 

o to analyze the results of the diffe- 
rent methods in light of the under- 
lying theories, thereby obtaining 

information on the most promising ways 

to achieve further improvement 

The 1980 Workshop had been successful 
in achieving these goals, so the new Work- 
shop was organized in a similar way. Two 
test cases were selected. The so-called 
HSVA tanker 2 was again chosen, as being the 
best documented test case available, see 
for instance Hoffmann 3 and Wieghardt and 
Kux 4 . Even though measurement data have 
been obtained only at model Reynolds num- 
bers, the participants were asked, as an 
opt ional exercise, to carry out calcula- 
tions also for a full scale Reynolds num- 
ber. This was to shed some light on the 
difficulties encountered at this larger 
scale, for which calculations have been 
very rare, so far. The second test case 
was designed to produce a significantly 
different stem flow with a minimum change 
of geometry , More information about the 
design philosophy is given in the next sec- 
tion. A very important feature of the se- 
cond case (the "Mystery case") was that no 
measurement data were available when the 
calculations were carried out. 

As in the first Workshop, attention 
was confined to double models, in which 
wave effects are absent and the free sur- 
face may be considered as a plane of sym- 
metry. Also, only the flow on the bare 
hull, without appendages and propulsors, 
was considered. 

The first announcement of the "1990 
SSPA-CTH- I IHR Workshop on Ship Viscous 
Flow" was distributed in late 1988, to- 
gether with a questionnaire to be returned 
by 15 May 1989 by all researchers interes- 
ted in taking part in the Workshop calcula- 
tions. Efforts were made to invite parti- 
cipation by originators of commercial CFD 
codes. No less than 21 computers indicated 
their Interest in participating, and in 
August the data and instructions for the 
first test case were sent out. Similar in- 
formation for the second case was distribu- 
ted in December. In the early summer of 
1990 results from 19 methods had been re- 
ceived. and the difficult task of conden- 
sing all the information into a uniform 
format was started at CTH. By the time of 
the Workshop meeting, which was held at CTH 
12-14 September 1990, all results had been 
plotted in a way such that comparison bet- 
ween the different methods could be easily 
made. Further, the replies to the question- 
naire had been compiled at IIHR and conden- 
sed into a table, useful for quick referen- 
ce to the theory behind each method. 

Zn the present paper a brief summary 
of the Workshop is given. The two test ca- 
ses are described next, followed by an 
overview of the methods. Thereafter, some 
important results are discussed and finally 
some conclusions are drawn. For a full ac- 


count of the Workshop, reference should be 
made to the Proceedings 5 . 


Test cases 
HSVA Tanker, Case I 

A body plan of the first test case, 
the HSVA tanker is shown in Fig. 1. The 
boundary layer measurements by Hoffmann 3 , 
and the subsequent stern- flow measurements 
by Wieghardt, Kux and Knaack 4 , were made on 
a double^ model of this hull in a 1.2 m dia- 
meter, slotted-wall wind tunnel, in which 
the turbulence level was of the order of 1 
percent. Different types of pressure probes 
were used* The model was supported in the 
tunnel by means of wires and a sting at the 
stern. The nominal length of the model was 
2.74 m but for reference length we have 
used the length between perpendiculars, L * 
2.664 m, whi ch gives a Reynolds number of 
5 • 106 . Neither the tunnel nor the support 
were modeled by any of the computers^ 

Mystery hull. Case II 

The second test case, for which no 
measured data were available by the time of 
the meeting, was designed by ProT G, Byne^ 
at Chalmers and Mr b.G. Jonsson at SSPA. 

The purpose was to create a significantly 
different wake pattern w ith a minimum of 
geometry change as compared to the first 
case. Thus, only the stern sections were 
modified. By making them more U-shaped 
stronger longitudinal vortices could be 
expected behind the hull, creating a more 
distorted wake field, see Dyne 6 . A body 
plan is shown in Fig. 2. 

Velocity measurements using Laser-Dop- 
pler Velocimetry were carried out at the 
University of Hamburg after the Workshop 7 , 
but the results are analyzed and included 
In the Proceedings 5 . 

Overview of methods 

Some 19 organizations from 12 count- 
ries participated in the Workshop. A sum- 
mary of the important characteristics of 
the methods is presented in Table 1 (at 
the end of the paper ) . This table was 
prepared on the basis of information sup- 
plied by the participants in a question- 
naire that was distributed at the beginning 
and again at the end of the Workshop. Ef- 
fort was made to obtain as much and as 
accurate information as possible on each 
method. The following is a review of some 
of the similarities and differences among 
the methods. 

The overall strategy summarized in 
item A indicates that 12 participants re- 
stricted their calculations to the stem 
and wake flow (S) while 7 treated the com- 
plete hull (H), including the bow. Both 
global (G) and zonal (Z) methods are repre- 
sented but the most common combination is a 
global method applied to the stem flow 
(S,G). The 5 zonal approaches employ an 
inviscid-flow method. In three of these. 




'on the other hand, satisfy the law of the 
wall and related conditions at one or more 
grid points away from the wall, and do not 
explicitly solve the flow equations in the 
near- wall layer. 

Section D of Table 1 provides an over- 
view of the boundary and initial conditions 
employed in the various calculations. As 
all methods approach the solution either by 
a time marching or an iterative process, 
all require initial conditions ( time - 0 ) . 
In this respect, three methods start the 
solution from rest, 7 use uniform flow, one 
uses boundary layer solutions, another uses 
a potential flow solution, while 5 indicate 
some other procedure. Among the last cate- 
gory are methods in which a parabolic march 
is made through the solution domain with an 
assumed initial pressure field. There is 
also considerable variation among methods 
in the quantities that are determined or 
prescribed at the start of the calcula- 
tions. The Initial conditions presumably 
Influence the number of iterations or time 
steps required to obtain the steady state 
solutions* that are sought. 

The upstream boundary conditions de- 
pend upon whether the solution is obtained 
for the entire hull or only the stem and 
wake flow. In the case of the calculations 
for the entire hull, two types of treatment 
have been made. In five of the seven such 
calculations, the upstream boundary is pla- 
ced about one-half ship length ahead of the 
bow, and quite simple boundary conditions 
are prescribed In the uniform flow there. 

In the remaining two (IL and KaO), the 
boundary layer over the bow is calculated. 
It appears that none of these calculations 
takes any account of the initial region of 
laminar flow or of the transition that was 
provoked artificially In the experiments. 

In the 12 calculations that were restricted 
to the stern and wake flow, there are dif- 
ferences in the location where the calcula- 
tions were started as well as in the para- 
meters that are prescribed. Three partici- 
pants started the calculations at a check 
station, X/L « 0.646, where the integral 
parameters of the boundary layer were 
supplied from experiment. (The X-axis is 
along the hull, with the origin at the 
bow.) Others started the calculations some- 
what ahead of this section, using. In some 
cases, the data at the check station as a 
guide. In all such calculations, however, 
the detailed velocity and turbulence para- 
meter distributions required by the methods 
had to be generated by the participants. In 
this regard, most appear to have used two- 
dimensional boundary layer correlations, 
with the three dimensionality neglected. 

In one case (PJL), however, a special up- 
dating scheme was devised to obtain a set 
of initial conditions that is consistent 
with the equations being solved. The diffe- 
rences in the initial conditions are likely 
to be observed most clearly in the results 
at the check station. 

There is also considerable variation 
in the location of the downstream boundary 


where the solutions terminate, and the con- 
ditions applied at that boundary. In most 
cases the solutions are taken far enough 
from the stern to assume a negligible up- 
stream influence and for the parabolic or 
extrapolated conditions to be valid. 

The boundary conditions at the hull 
surface were discussed above in connection 
with the turbulence model. For complete- 
ness, however, we note that the no-slip 
conditions are applied explicitly in some 
methods whereas they are satisfied indi- 
rectly in methods that rely on the wall- 
f unctions approach. 

All of the calculations presented at 
the Workshop have exploited the geometric 
symmetry about the ship centerplane and 
calculated only one half of the hull. Also, 
all participants assumed a double body and 
applied symmetry conditions along the water 
plane. 

The location of the solution boundary 
In the "farfield" some distance from the - 
ship axis, and the conditions specified 
along that boundary, also show considerable 
differences. With respect to the location, 
the five zonal calculations performed with 
viscous and inviscid methods use a boundary 
that varies from 0.08 of a ship length in 
one case (H) to 0.7 in another (KaO). Re- 
call that this is the boundary at which the 
two solutions are matched and, as noted 
earlier, none of the calculations accounted 
for the interactions between the two 
regions. In these zonal calculations, the 
inviscid solutions provided the boundary 
conditions for the viscous solutions. In : 
the remaining, global approaches, the loca- 
tion of the farfield boundary ranged from 
0.1 of a ship length (T) to 1.5 (HC) with 
many using a value of one ship length. 
However, those who have used boundaries 
rather close to the hull (SZC and T) have 
provided boundary conditions from Inviscid 
flow. In this respect, these methods, cha- 
racterised as global by their originators, 
could also be regarded as zonal. Methods 
that have placed the boundary at larger 
distances from the hull have tended to pre- 
scribe uniform-flow conditions In the far 
field. Some differences are observed, how- 
ever, in these cases in the particular 
variables or conditions that are specified 
or satisfied; see, for example, codes C and 
A in the table. 

Section E of Table 1 pertains to the 
generation and control of the numerical 
grid. For the Workshop calculations, all 
participants employed a single block grid 
although one (GCHM) Indicated that their 
method can accommodate a multi-block grid. 
Fully nonorthogonal grids as well as grids 
that are orthogonal In some planes (typi- 
cally In the transverse sections) or at 
boundaries (usually at the hull ) are rep- 
resented in the calculations. The most po- 
pular method for generating the grid appe- 
ars to be numerical, although some partici- 
pants have employed analytic and algebraic 
methods, and even combinations of methods. 



Among the numerical approaches, elliptic 
methods for the entire three-dimensional 
solution domain are the most common. Seve- 
ral methods use post-generation smoothing 
of the grid. Control of the grid is exer- 
cised most commonly from the boundaries 
although some indicate that it is done from 
inside the solution domain. 

The numerical features of the various 
methods are summarized in section F. First 
of all, we note that there is no participa- 
tion from the finite-element community and 
therefore items in the questionnaire that 
were designed specifically to obtain infor- 
mation on such methods are deleted in the 
table. Of the methods represented at the 
Workshop, there are four finite-difference 
methods, 10 finite- volume methods, and, for 
lack of a better term, four are classified 
as mixed methods. To these must be added 
the one based on integral equations which 
are also solved by finite-difference 
methods. It should be pointed out that not 
all of the methods in each of these catego- 
ries are generically different. In other 
words, there are groups of methods in each 
class which share a great deal in common 
and, from a numerical perspective, may be 
classified as a single method. Be that as 
it may, we shall note the most significant 
features of the methods represented here. 

First, discretization in both stag- 
gered and regular (colocated) grids is used 
in finite-difference as wej.1 as finite- 
volume methods. There is some correlation 
between the type of grid used and the met- 
hod employed to solve for the pressure (or 
establish the pressure-velocity coupling) 
for this incompressible flow. Differencing 
of the convective terms is made by a vari- 
ety of means, including upwinding, central 
differences, hybrid combinations and use of 
analytic solutions (FA). The responses to 
the questions on formal order of accuracy 
and formally conserved quantities are 
rather surprising. Most methods claim first 
or second-order accuracy but two consider 
their methods to be accurate upto the third 
order. The response to the second question 
seems to be correlated with the equations 
that are solved rather than any formal 
attempt in the methods to conserve mass, 
momentum, and total energy. 

For the pressure, four methods use 
fully-coupled solutions. Two of these (GCHM 
and K) employ artificial compressibility, 
one (AS) uses a discretized continuity 
equation, another (H) uses the normal mo- 
mentum equation. Most of the remaining 
methods employ segregated pressure-velocity 
coupling algorithms, SIMPLE being the most 
common • 

In the execution of the solutions, 
most methods employ iterative techniques 
with under-relaxation or variable time 
steps. Explicit as well as implicit methods 
are represented. The solutions of the dis- 
cretized algebraic equations are obtained 
using different tactics, including point 
substitution, line substitution, LU decom- 
position of matrices, and ADI methods. 


Section G of Table 1 concerns the com- 
putations performed by the participants. As 
the two test cases are not substantially 
different with respect to the computational 
effort Involved, the numbers In this sec- 
tion are typical of both cases. First of 
all, it is quite significant that calcula- 
tions have been performed not only on 
state-of-the-art supercomputers (designated 
by S) but also on smaller machines (desig- 
nated by M), such as workstations. The to- 
tal number of grid points employed shows a 
great deal of variation, ranging from a low 
of 8,000 to a high of 253,000. It is Inte- 
resting to note that the higher numbers are 
not necessariliy correlated with the use of 
supercomputers, nor are they correlated 
with the use of near-wall turbulence models 
or calculations made for the entire hull 
including the bow. 

The number of time steps or iterations 
performed to obtain the solutions presented 
here varied widely, ranging from only 40 
(BZLS and SZC) to 20,000 (ZM), although, in 
the latter case, the large-eddy simulation 
should never realize a steady state as de- 
fined for the other methods. It should be 
noted that a variety of convergence cri- 
teria were adopted to declare that a steady 
state had been obtained. The storage requi- 
rements and computer run times also varied 
over very wide limits. This is not surpri- 
sing in view of the differences in the num- 
ber of grid points and machines employed, 
but it is interesting to note that run 
times ranging from several hours to 5 days 
were reported by users of the smaller ma- 
chines. The differences in the machines are 
also reflected in the cpu time per itera- 
tion per grid point. The fastest times were 
of course reported by users of supercompu- 
ters. 


Results and Discussion 

Results were requested from the parti- 
cipants at a Reynolds number of 5 • 10 6 for 
both hulls. As an optional exercise the 
computors were also asked to submit full 
scale results corresponding to a Reynolds 
number of 2 • 10 9 for case I. All but one 
delivered the model scale predictions but 
only three had computed the high Reynolds 
number case. 

Pressure and friction distributions 
were reported at the waterline and keel and 
along three section girths on the hull, 
while velocity distributions in the form of 
iso- velocity contours and, in some cases, 
cross- flow vectors were given at four sec- 
tions. For the most Interesting section, 
the propeller plane, the pressure and (if 
computed) the turbulent kinetic energy were 
also reported. 

The complete results may be found in 
the Proceedings 5 , but in the present paper 
only a few representative examples will be 
given. These include the pressure distri- 
bution along three section girths and the 
velocity distribution (iso-velocities and 
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cross-flow vectors) at the propeller plane 
of Case I , see Figs. 3, 4 and 5. In these 
Figures the results of all 19 methods are 
included. Measured data are represented by 
symbols in Fig. 3, while the computed re- 
sults are given as lines. In Fig. 4 iso- ve- 
locity contours corresponding to 0.3, 0.4, 
0.5, 0.6, 0.7, 0.8 and 0.9 times the undis- 
turbed velocity are given. For comparison 
the experimental data are shown in the lo- 
wer right corner. Note the different scales 
for the computed and experimental cross - 
flow vectors in Fig. 5. 

Before embarking on a more detailed 
analysis a comparison of the results pre- 
sented at this Workshop with those obtained 
a decade ago for the SSPA-ITTC Workshop on 
Ship Boundary Layers 1 is quite revealing. 
Fig. 6 shows the axial velocity contours 
for the HSVA Tanker at X/L-0.942 presented 
at the previous Workshop. A comparison of 
these with the results shown in Fig. 4 pro- 
vides an overview of the achievements of 
the past decade. It is clear that the ear- 
lier boundary layer methods have given way 
to those based on the Navier-Stokes equa- 
tions. Only one such method was represented 
in 1980. At the present Workshop, only one 
boundary layer method was represented. The 
question is whether or not real progress 
has been achieved in the prediction of the 
flow. If the contours of axial velocity at 
the stem are used as the only measure of 
success, then we may conclude that progress 
has indeed been made. But, consider the 
following observation. Most calculations 
methods of the past did rather well at pre- 
diction of the boundary layer over the hull 
and failed only at the stern. Among the 
present methods, few, if any, predict the 
boundary layers with the same level of ac- 
curacy but continue to provide results for 
the flow over the stem and into the wake. 

The girthwise pressure distributions 
of Fig. 3, will be considered section by 
section. First of all, we note that there 
are significant differences in the way dif- 
ferent methods obtained the results at the 
check station, X/L - 0.646. As indicated in 
item A of Table 1 (Code H), seven methods 
performed calculations for the entire hull, 
starting upstream of it. Although the man- 
ner in which this was accomplished dif- 
fered, their results at the check station 
reflect, to some extent, their resolution 
of the flow over the bow. Of the remaining 
methods, three started the calculations at 
the check station itself (see item D, Table 
1) using the integral parameter information 
provided there. Note that no information on 
the crossflow at the check station was pro- 
vided. The remaining methods started the 
solutions on the hull somewhere upstream of 
the check station and may have used the 
information provided to guide the selection 
of the upstream conditions. These differen- 
ces among the methods must be taken into 
consideration when examining the results at 
the check station as well as further down- 
stream. 

In view of the foregoing, the results 
of four of the seven methods that calcula- 



6. Axial velocity contours at X/L * 0.942, 
predicted in the 1980 Workshop 1 . 

ted the flow over the entire hull are in 
remarkably good agreement with the data 
with respect to the girthwise pressure dis- 
tribution at the first station. In fact, 
their predictions are as good as those of 
some of the methods that started the solu- 
tions on the hull. Methods that started the 
calculations on the hull show varying le- 
vels of agreement in the pressure distribu- 
tion at the check station. The disagreement 
at the check station of methods that star- 
ted at that station, or just upstream, are 
rather surprising. The reasons for this are 
not clear. Be that as it may, these diffe- 
rences should be borne in mind as compari- 
sons are made further downstream. 

At the next downstream section, X/L 
- 0.875, we see that the measured pressure 
distribution shows a decrease from the 



keel, with a minimum Cp of about -0.24 
around 30 percent of tne girth, an increase 
up to about 70 percent of girth, followed 
by a near constant value of -0.11 around 
the waterline. It is clear from Fig. 3 that 
many methods reproduce this trend but, with 
one possible exception, fail to predict 
either the magnitude or the location of the 
pressure minimum. 

At the last section, X/L - 0.942, the 
measurements indicate minimum cp at around 
8 percent girth, much closer to the keel 
than at the previous station, and an in- 
crease followed by a constant pressure 
around the waterline. The calculations ge- 
nerally predict higher pressures throughout 
and all fail to capture the dip in pressure 
around 45 percent girth. The_ methods that 
did poorly at the previous station continue 
to yield poor results at this station also. 

From the point of view of propeller 
design the most interesting results are 
those at the propeller plane, reported in 
Figs. 4 and 5. The propeller disk is indi- 
cated in these Figures, and it is seen that 
particularly in this area the results vary 
considerably between the different methods. 
For classification purposes it is possible 
to distinguish between three types of re- 
sults, see Fig. 7. 


ABC 


7. Different types of axial velocity con- 
tours at the propeller plane. 

A. V-shaped contours 

B. Contours with a bulge at or below the 
propeller center, indicating that part 
of the ship boundary layer has been 
displaced downwards by the longitudi- 
nal vortex hitting the propeller plane 
slightly above the propeller center 
(see Fig. 5). 

C. Contours with a pronounced "hook" due 
to the action of a stronger vortex. 

Obviously the measured wake is of type 
C with quite distinct hooks in the 0.3 and 
0.4 contour curves. 

Investigating the contours of Fig. 5 
it is seen that only one method produced a 
wake of type C, but the results of this 
method outside the propeller disk, particu- 
larly near the waterline, look quite unre- 


alistic. Four or five methods predicted a 
wake of type B, while the remaining ones 
produced A type wakes. 

For a propeller designer this situa- 
tion is unsatisfactory, since the wake dis- 
tribution determines the variation in loa- 
dina durinq one turn of the propeller, i.e. 
the vibratory forces. On the other hand, as 
was shown by one of the authors (Larsson) 
at a continued Workshop in Osaka in the 
fall of 1991, the mean value around the 
circle at all radii may still be well pre- 
dicted, as may the mean value of the velo- 
city in the entire propeller disk. This 
means that the pitch and its distribution 
on the blades might be reasonably well pre- 
dicted using the calculated wakes of the 4- 
5 best methods. 

The reason for the failure to predict 
hook-like velocity contours was discussed 
extensively during the Workshop. To struc- 
ture the discussion it was divided into 
four main themes: grid generation, equa- 
tions and boundary conditions, turbulence 
models and numerical method. 

There appeared to be a general con- 
sensus that grid resolution was not a ma- 
jor cause of the differences between, on 
the one hand the methods themselves, and on 
the other hand the measurements. This opi- 
nion was substantiated by the fact that 
several computors had carried out grid in^ 
dependence studies with very small changes 
in the results. There was however a general 
feeling that the resolution of the trans- 
verse pressure gradient was too low near 
the bilge (region of high transverse curva- 
ture) on the afterbody. Since this is where 
the longitudinal vortex is generated, the 
lack of resolution could explain the too- 
weak vorticity in the wake predicted by 
most methods. An impression of the grids 
used can be obtained from Fig 5, although, 
for clarity, not all the grid points are 
shown for some methods. 

As for the governing equations, a dif- 
ference not reported in the questionnaire 
turned out to be the way in which the tur- 
bulence terms are considered. It became 
clear that only a few methods include all 
of the terras. Another difference in the 
governing equations (appearing from Table 
I) is that some methods use the partially- 
parabolic approximation, while others are 
fully elliptic. This difference was discus- 
sed, and It was concluded that no correla- 
tion could be found between the approxi- 
mation level in this respect and the per- 
formance of the method, as judged by the 
results reported. Some participants had in 
fact computed Case I using both types of 
equations and found very small differences. 

Different inlet boundary conditions 
are required depending on whether the com- 
putational domain starts on the hull or 
upstream of It, but in all cases the parti- 
cipants were free to match their solution 
to the measured data at X/L - 0.646. Large 
differences are, however, seen at this sta- 
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tion, indicating that this possibility was 
not exploited by most of the computors. 
Instead, the boundary layer at the inlet 
station, if on the hull, seems to have been 
estimated from flat-plate correlations. The 
methods starting upstream use uniform inlet 
flow. To avoid numerical problems a nonzero 
value of the eddy viscosity had to be app- 
lied right from the start. Values of 50 to 
100 times the laminar viscosity were men- 
tioned. No attempt was made to consider 
transition. The general feeling was that 
the inlet conditions did not significantly 
influence the results at the stations on 
the afterbody and in the wake. This conclu- 
sion had in fact been verified by some par- 
ticipants . 

Some discussers expressed the opinion 
that the only way to resolve details of the 
flow close to the wall ( such as limiting 
streamlines) is to abandon the wall law. 

The general opinion was, however, that it 
is difficult to see a large difference in 
performance between the methods using the 
wall law and the others. A better predic- 
tion of the wake contours would have been 
expected, and this was achieved by some of 
the non-wall law methods but not all. It 
was pointed out also that virtually the 
only way to compute the full scale case is 
to employ the wall law, since otherwise the 
innermost grid points have to be positioned 
extremely close to the surface, giving rise 
to problems with the cell aspect ratios 
(the non-dimensional distance from the sur- 
face, y/L, is about 250 times smaller for 
the full scale case for a given value of 
y+ ) . The other boundary conditions were 
considered relatively unimportant for the 
problems at hand. 

Different opinions on the general im- 
portance of turbulence modeling were ex- 
pressed. It was argued that as experiments 
have Indicated that the Reynolds stresses 
are very small in the major part of the 
viscous region near the stern, even an in- 
vlscld calculation might produce a reaso- 
nable result. Some participants reported on 
earlier computations for axisymmetric and 
three-dimensional bodies, where this app- 
roach had been tested with relatively good 
results. Obviously the inviscid region had 
to be restricted to the neighborhood of the 
stem. 

A possible explanation for the failure 
to predict the correct wake contours might 
be the inability of the methods to resolve 
the pressure field accurately, i.e. to pre- 
dict the transverse pressure gradients that 
are believed to produce the vortical flow 
structure. In this case, a number of nume- 
rical Issues are involved, including the 
grid arrangement used (staggered versus 
regular) and the manner in which the pres- 
sure is calculated* 


Conclusions 

The Workshop clearly showed that great 
progress has been made during the 1980 T s in 
the development of methods for predicting 


the flow around the stern and in the near 
wake of ships. The boundary layer based 
approaches of the 1980 Workshop have given 
way to methods based on the Reynolds -Ave- 
raged Navier-Stokes equations, albeit using 
relatively simple turbulence models. While 
the former methods either broke down be- 
fore, or predicted completely unrealistic 
results in the propeller plane, the new 
methods in general capture the gross featu- 
res of the wake, such as the thin shear 
layer in the lower part and the piling- up 
of boundary layer low speed flow around 
half draught. In fact, the best methods of 
the Workshop predict the shape and location 
of the velocity contours in most of the 
propeller plane with reasonable accuracy. 
The results are however less satisfactory 
in the central part of the wake, i.e. in- 
side the propeller disk. The bilge vortex, 
although weak, redistributes the low speed 
flow from the boundary layer in such a way 
that very uneven hook shaped velocity con- 
tours are created. This feature is missed 
to varying degrees by the different 
methods. Various reasons for this were dis- 
cussed during the Workshop, but no definite 
answer could be given. One possible expla- 
nation is that the large velocity and pres- 
sure gradients in the bilge region are too 
inaccurately resolved in the grids employ- 
ed. 

There are no general differences in 
performance between methods based on zero 
equation turbulence models as compared to 
the two equation models. The disadvantages 
of the simpler models may be outweighed by 
the advantages of computing the flow all 
the way down to the surface. Two equation 
models are usually employed in combination 
with the wall law. 

Although the new methods are superior 
to the old ones in predicting the wake, 
results from calculations starting upstream 
of the hull, or in front of the check 
station. Indicate that the ability to pre- 
dict the thin boundary layer has not been 
improved, rather the contrary. The likely 
reason for this is grid resolution. To save 
computer time very few grid points are lo- 
cated in the thin boundary layer on the 
forebody, where the much faster boundary 
layer methods may use grids with a very 
high resolution. This suggests a zonal app- 
roach where the expensive Navier-Stokes 
method is used in the stem and wake region 
while an efficient boundary layer method is 
used for the rest of the hull. 

A final point to note is that most 
methods predicted, at least qualitatively, 
the differences between the two test ca- 
ses. The computed cross flow for the second 
case is considerably larger in the propel- 
ler disk than for the first case, as the 
measurements have indicated, and the change 
in the contours of axial velocity also 
shows the right trend. 
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Segregated: Poisson equation, e.g., capacitance techniques (P); pressure correction 
technique, e.g., SIMPLE, SIMPLER, SIMPLEC. PISO, saddle-point method, 
preconditioning (C, list); other (e.g., boundary -layer type code) describe (O) 

Code - S,; SIMPLE; S* SIMPLER; S* SIMPLEST; S 4 : SIMPLEC; P: PISO 







